NAA 
IT] 


Vol. 39 


HIT 


o. 8 


第 39 AR 8 期 计算 机 应 用 研究 
录用 定稿 Application Research of Computers Accepted Paper 
pp p 


ETZERA Z d dE ZAETERETS EHE C 


» NIE. y T 
4b oo, prj o 


(1. 东南 大 学 自动 化 学 院 , 南京 210096; 2. 复杂 工程 系统 测量 与 控制 教育 部 重点 实验 室 ， 南 京 210096; 3. 北京 航空 
航天 大 学 杭州 创新 研究 院 ， 杭州 310051; 4. 北京 航空 航天 大 学 ,北京 100191) 


H 要 : 针对 多 入 多 出 (Multiple Input Multiple Output, MIMO) 非 线 性 时 滞 系 统 辨 识 的 准确 性 与 实时 性 问题 ， 提 出 
基于 多 维 泰勒 网 (Multi-dimensional Taylor Network, MTN) 的 辨识 方案 。MTN 作为 辨识 模型 ， 综 合 利用 权 剪 枝 
(Weight-Elimination, WE) 算法 和 共 罗 梯度 算法 (Conjugate Gradient, CG), PP WE-CG 算法 作为 MTN 办 识 模型 的 学 
习 算 法 ; WE 算法 可 以 有 效 精简 MTN 辨识 模型 结构 ， 从 而 降低 计算 复杂 度 ， 提 高 模型 实时 性 能 。 最 后 ， 引 入 一 个 数 
值 仿真 例子 和 一 个 工程 实例 来 验证 所 提 辨 识 方案 的 有 效 性 ， 同 时 与 传统 的 MTN 辨识 方案 作对 比 ， 给 出 了 二 者 的 准 
确 性 与 复杂 度 分 析 ， 突 出 所 提 闪 识 方案 的 准确 性 与 实时 性 。 实 验 结果 表明 ， 所 提 方案 够 准确 地 对 MIMO 非 线 性 时 灌 
系统 进行 辨识 。 同 时 ， 相 比 传统 的 MTN 办 识 方案 ， 所 提 辨 识 方案 结构 更 精简 ， 具 有 更 低 的 算法 复杂 度 。 

关键 词 : 多 维 泰勒 网 ; 办 识 模型 ;多 入 多 出 系统 ; 非 线 性 时 灌 系 统 ; 剪 枝 算法 

中 图 分 类 号 : TP271 doi: 10.19734/j.issn.1001-3695.2022.01.0031 


Identification of MIMO nonlinear time-delay system using multi-dimensional Taylor network 
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Complex Systems of Engineering, Ministry of Education, Southeast University, Nanjing 210096, China; 3. Hangzhou 
Innovation Institute, Beihang University, Hangzhou 310051, China; 4. Beihang University, Beijing 100191, China) 


Abstract: For the accuracy and real-time performance of multiple input multiple output (MIMO) nonlinear time-delay system, 
this paper propose the identification scheme based on multi-dimensional Taylor Network (MTN) . It use MTN as an 
identification model. MTN identification model's learning algorithm adopts the WE-CG algorithm, which combines the 
weight-elimination (WE) algorithm with the conjugate gradient (CG) algorithm. WE algorithm can simplify the structure of 
MTN model so as to reduce computational complexity and improve the real-time performance. Finally, there are two 
experimental examples containing a numerical example and a project one to verify the effectiveness of the proposed 
identification scheme. And use the traditional MTN identification scheme to compare. Both the accuracy and the 
computational complexity analysis illustrate the accuracy and the real-time performance of the proposed scheme. Results 
from the experiments and comparison show that the proposed scheme can identify the MIMO nonlinear time-delay system 
accurately, and the proposed identification scheme has a simpler structure and lower computational complexity than the 
traditional MTN identification scheme. 

Key words: multi-dimensional Taylor network; identification model; MIMO system; nonlinear time-delay system; weight- 
elimination algorithm 


神经 网 络 的 辨识 方案 ; 文献 [9] 针 对 非 线 性 系统 辨识 ， 提 出 了 
一 种 新 型 神经 网 络 方案 ; 文献 [10] 针 对 非 线性 动态 系统 , 提出 


0 引言 


工业 系统 中 大 都 具有 非 线 性 、 时 灌 等 特性 ， 给 建 模 和 控 了 一 种 基于 数据 知识 的 模糊 神经 网 络 辨识 方案 。 然 而 ， 
制 带 来 挑战 0。 以 往 在 线性 系统 中 ， 因 为 可 以 得 到 其 精确 的 。 神经 网 络 由 于 在 训练 过 程 中 花费 时 间 长 ， 易 陷入 极 小 值 
数学 模型 ， 所 以 利用 传统 的 辨识 方法 可 以 满足 其 要 求 ， 比 如 。 点 ， 且 算法 复杂 度 较 高 ， 从 而 不 能 满足 实时 构建 预测 模 
最 小 二 乘法 、 阶 跃 响应 法 等 B 轨 。 而 对 于 非 线性 时 滞 系 统 辨识 “型 的 要 求 00。 此 外 ,， 相 比 于 单 入 单 出 (Single Input Single 
问题 ， 由 于 很 难 获得 系统 精确 的 数学 模型 ， 因 此 ， 以 往 针 对 “Output，SISO) 系 统 ， 多 入 多 出 (Multiple Input Multiple 
线性 系统 的 辨识 方法 很 难 达 到 理想 的 效果 ， 从 而 给 辨识 带 Output, MIMO) 系统 由 于 输入 和 输出 变量 之 间 的 关系 
困难 所 。 随 着 神经 网 络 的 出 现 ， 由 于 其 强大 的 非 线性 逼近 能 ”更 加 复杂 而 更 加 难以 建 模 02 -14 。 


H 


力 ， 给 非 线性 辨识 提供 了 一 种 有 效 的 解决 途径 [9I0。 比 如 ， 针对 以 上 问题 ， 严 洪森 教授 于 2010 年 10 月 提出 多 维 泰勒 
文献 [6] 针 对 非 线性 系统 ， 提 出 基于 人 工 神经 网 络 与 切 比 雪夫 网 的 建 模 思 想 0549， 之 后 被 应 用 在 非 线性 时 间 序 列 建 模 与 故障 
函 数 的 时 滞 自 适应 辨识 方案 ; 文献 [7] 基 于 时 滞 神 经 网 络 与 核 监测 上 07-20。 之 后 在 2010 年 12 月 严 洪森 教授 提出 多 维 泰勒 网 
函数 的 内 部 关系 ,提出 针对 工程 领域 的 非 线 性 系统 辨识 方案 ; 优化 控制 的 思想 0549， 应 用 于 非 线 性 系统 控制 领域 P2 约 ， 先 后 
文献 [8] 针 对 带 有 未 知 时 灌 的 非 线 性 连续 系统 ,提出 基于 时 灌 “对 非 线性 时 变 系统 234、 非 线性 时 滞 系 统 P3] 与 非 线性 随机 系统 24 
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上 


的 控制 问题 进行 了 广泛 研究 。 多 维 泰勒 网 与 多 维 泰勒 网 优化 
控制 器 县 有 以 下 特点 05 09:， a) MTN 具有 可 以 任意 逼近 非 线 
性 函数 的 能 力 ; b) MTN 模型 与 控制 器 ， 其 结构 简单 ， 只 包含 
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函数 任意 准确 地 逼近 。 
引 理 203 对 于 定义 于 一 个 闭 区 间 的 连续 函数 f(%,%,…,*) ， 


加 法 和 乘法 ， 不 含有 超越 函数 和 无 理 函数 ， 参 数 调 整 方便 ， 
便于 工程 实现 ， 无 须 复杂 的 数学 计算 ， 便 于 实现 实时 建 模 与 
控制 ，c) MTN 控制 器 本 身 具有 “学 习 ” 能 力 ， 使 它 在 控制 过 
程 中 不 断 地 完善 自己 ， 使 控制 效果 越 来 越 好 ， 提 高 控制 器 的 


En 


RT DLE Ep 逼近 。 其 中 Nm 为 逼近 展开 式 中 乘积 项 的 总 


项 数 ， 风 为 逼近 展开 式 中 第 :个 乘积 项 之 前 的 权 值 ，4, 为 展 
开 式 中 第 :个 乘积 项 中 变量 * 的 肾 次 。 


正如 图 1 所 示 , MTN 采用 前 向 单 中 间 层 结构 , 包括 输入 


鲁 棒 性 能 ; d) MTN 控制 器 具有 良好 的 可 扩展 性 , 能 够 融合 其 
他 控制 方法 或 者 控制 思想 ， 具 有 一 定 的 容错 性 和 进化 能 力 。 

在 基于 此 ， 本 文 针对 MIMO 非 线性 时 滞 系 统 ， 提 出 基于 
MTN 的 辨识 方案 。 首 先 ，MTN 作为 辨识 模型 ， 然 后 ， 综 合 
F FL BE d (Weight-Elimination, WE) £t YE JU 44 90 1% BE $2 1: 
(Conjugate Gradient, CG)， 即 WE-CG 算法 作为 MTN 辨识 模 
型 的 学 习 算 法 ，WE 算法 对 MTN 辨识 模型 结构 进行 精简 从 
而 降低 计算 复杂 度 以 满足 实时 性 要 求 ， 最 后 ， 引 入 一 个 数值 
仿真 例子 和 一 个 工程 实例 -连续 搅拌 反映 器 (continuous stirred 
tank reactor, CSTR) 来 验证 本 文 所 提 方 案 的 有 效 性 ， 同 时 引入 
传统 的 MTN 辨识 方案 07 作 对 比 ， 给 出 了 二 者 的 准确 性 与 复 
杂 度 分 析 。 

本 文 的 主要 贡献 如 下 ; 

a) 针对 MIMO 非 线性 时 滞 系 统 , 借助 MTN 的 强大 逼近 
性 能 ， 提 出 基于 MTN 的 辨识 方案 ， 所 提 方 案 易 于 实施 ; 
b) WE-CG 算法 作为 MTN 辨识 模型 的 学 习 算 法 ，WE SE 
法 可 以 有 效 精 简 MTN 辨识 模型 结构 ,从 而 降低 算法 复杂 度 ， 
提高 模型 实时 性 能 ; 
c) 所 提 方 案 具 有 较 低 的 计算 复杂 度 ， 即 实时 性 能 较 好 ; 
d) 所 提 辨 识 方 案 的 具有 一 般 性 和 实际 应 用 价值 。 


1 ”系统 描述 与 多 维 泰勒 网 


1.1 系统 描述 
MIMO 离散 非 线 性 时 滞 系 统 由 如 下 表示 : 
Ik) - gly(k 2). y(k 2), y(k n,), 0 
u(k—d),---,u(k-n, —d)] 
其 中 ， gO-LaO.-s e, OI eR" 为 未 知 的 非 线 性 向 量 函 数 ; 
YW 2 [y, (9... y, QI ER; uk) =u K), ER 分 别 为 系统 (1) 
的 a 维 输出 和 P 维 输入 ; mw ，(j=12.…q )，m (i=12…:p) 分 
别 为 第 j 个 输出 分 量 和 第 i 个 输出 分 量 的 阶 次 ， 目 


可 


层 、 中 间 层 和 输出 层 。 同 时 基于 引 理 1 和 引 理 2 说 明 ， 只 要 
Nm 足够 大 ， 它 就 能 够 以 足够 的 精度 逼近 任意 模型 0825。 


gk) rh (k+l) 
£0) » X) 
&( EX 

», +1) 


N $ p AN 
Neap d 


UU X : 


图 1 MIMO 的 MTN 结构 图 

Fig. 1 MIMO MTN structure 
1.2.3 多 维 泰 勒 网 复杂 度 分 析 
于 使 用 相同 的 计算 机 ， 相 同 的 时 间 来 做 实时 性 的 分 析 
是 不 可 能 的 ， 因 此 本 文 引用 浮 点 数 来 估计 方案 的 计算 量 ， 即 
实时 性 2627。 一 个 加 号 和 一 个 乘 号 都 叫做 一 个 浮 点 ， 加 号 和 
乘 号 的 总 数 作为 总 的 浮 点 数 ， 即 浮 点 数 的 总 和 代表 着 计 
算 复杂 度 。 当 展开 次 军 为 2 次 ， 表 1 MTN d — Vs 
JP B VE SER AS HEUT, RH MTN 因 其 结构 简单 而 带 来 
出 色 的 实时 性 能 。 


n, = (n,,n,,,...,n,.] 


, n.—Íngnog224n4) . 


d 为 系统 的 时 灌 。 
则 ， 在 上 上 时刻， 系统 (1) 的 第 j 个 输出 分 量 可 以 表示 为 
yik) 2 2jOn(k - D. y (e =n, ) 
Ya lk D, y, (kn, ), 
uj (k-d),u(k d —V,--, (2) 
u(k —d —n, ),:,u,(k-d), 
u,(k -d D, ,u,(k-d -n,) 


1.2 多维 泰 勒 网 
1.2.1 MIMO 多 维 泰 勒 网 

针对 于 MIMO 非 线 性 系统 ， 基 于 多 维 泰勒 网 ， 第 7 个 输 
出 分 量 多 +D 可 写成 如 下 形式 (13 


AG an) NM 

$,k«no 2 Ww do p 8? o0 Q) 
其 中 ， NEm RR P RA AIO JC X i EEE RRR 
Hs XC, CER F8] ERE 9; 00 m I9 00e uu COT, 100 表示 式 
(3) 中 第 7 个 乘积 项 的 权重 系数 ，h,5) 表示 第 ?个 乘积 项 中 变 

E 25Co 的 次 数 ， 且 OE : 
MIMO 非 线 性 系统 的 MTN 结构 医 
1.2.2 多 维 泰 勒 网 逼近 性 能 


IS 
B 


图 1 所 示 。 


表 1 MTN 的 计算 复杂 度 


Tab. 1 Computation complexity of MTN model 
操作 项 MTN 
中 间 节 点 数 n(n —/2-- 2n, +1 
TUB n, Oi (n D/2--2n, +1) 
加 法 n no no —1)/2 +27) 
乘法 n, (m(n —1)/2+2m +1) 
fO n 


RP n Mn SARR MTN 输入 和 输出 层 的 节点 数 。 
MTN 只 包括 加 法 与 乘法 , 其 计算 复杂 度 仅 相当 于 神经 网 络 中 
单个 神经 元 的 泰勒 展开 。 


2 MIMO 非 线性 时 滞 系 统 多 维 泰勒 网 辨识 方案 


本 文采 用 MTN 作为 辨识 模型 ，WE-CG 算法 作为 学 习 算 
ik, WE 算法 可 以 有 效 精简 MTN 网 络 结构 。 
2.1 MTN 辨识 模型 


对 于 系统 (1)，sg[] 可 以 由 MTN 以 任意 精度 逼近 ， 其 映射 
关系 记 为 801, ， 则 MTN 辨识 模型 可 以 写成 如 下 形式 : 
$00 2 &Dy(k - D. y(k 22). y(k-n,), 7 


u(k — d), --,u(k ^ n, ^ d), v4(k)] 


引 理 1051 任何 定义 于 一 个 闭 区 间 的 连续 函数 可 以 用 多 项 式 


Hip, $00 Jy MTN 辨识 模型 输出 , EL 900 9D 00. $000 ER ; 
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900 为 MTN 权 系 数 向 量 ; 6O - O 
m, n d 与 式 (1) 有 相同 的 含义 。 


然后 ， 系 统 (4) 在 第 xt 时刻 的 第 /次 输出 分 量 为 L 


&,Ol' eR", 
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为 了 获得 一 个 精确 的 模型 ， 
过 数据 fwelo,yK+D} 来 更 新 。 梯 度 法 是 一 
可 以 对 参数 进行 调整 ， 即 ， 


ulk) , y) ， 


权重 系数 向 量 疡 的 


», 可 以 根据 每 次 离线 
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值 应 该 通 


种 常用 的 优化 方法 ， 
MTN 的 


3,00 = 8L 7D, y on, ) s y TD ff ER SEZT RIDE HEURE K. B. Jy o A c 步 训 练 后 的 值 ， 则 
y,(k n, ),u(k - d), u(k d -V,---, E . m E 8$. P — 
Lid sa ada i tcd edis () Pash -AE , IH OT ue o AM, 梯度 法 采用 锯齿 形 的 
uk -d n, wj] 方式 到 达 最 小 值 点 ， 所 以 其 搜索 方向 会 一 直 保 持 垂 直 。 幸 运 
为 了 方便 而 不 失 一 般 性 ， 记 的 是 ， 采 用 共 斩 梯 度 法 可 以 有 效 地 解决 这 一 问题 ， 其 权重 可 
Loy uad 以 更 新 如 下 : 
可 Wan = wr +A P. (11) 
M n "n _ 81g. ^ ê: 10p, 
t= ( y*D H , Q——EStSaDQa , dT; "E ^A " 5 
Ye fix Bre deo pk "Uy DRUG 
[ETE B --£& 
&0O 13,00. —-8, 00.5, (0. 500]! = dic E NEM » 
" 随 着 MTN 的 输入 节点 增加 , 会 导致 “ 维 数 爆炸 ”， 因 此 
[y -D, y, (k =n, ) un (k-0), us (c-r, —d)] Piece dul Rb d 
有 必要 精简 网 络 结构 。 权 剪 枝 算 法 可 用 于 精简 MTN 的 中 间 
MTN 的 第 j 次 输出 分 量 I 的 权 值 向 量 为 节点 , 以 此 来 减 小 网 络 的 复杂 度 。 本 节 通 过 结合 CG 法 与 WE 
9,09 2E GO... is OT ik, BU WE-CG 法 作为 MTN 辨识 模型 的 学 习 算 法 。 
则 ， 相 应 的 MTN 模型 (5) 可 以 改写 成 如 EU 引 理 309] 假设 目标 函数 wo 可 由 最 小 结构 网 络 Net 
AG n É x 实现 ， 网 络 Net, 比 Netan 多 一 个 中 间 节 点 。 ER. Net, 能 实现 
$,00- 3 e, o [2/00 (b 7A ee dy d ' 多 i bd UET MERA 
f p h(x), 则 Net, 或 者 j h-1 个 中 [R] "HW 点 与 中 [R] WH 点 Neto 元 
Hr, Nun) 表示 ?元 函数 甸 0) 展开 成 交 次 通 近 多 项 式 乘积 项 EA, 而 余下 的 那个 独立 中 间 节 点 输出 权 值 为 零 ; 或 者 
WADA, WO 表示 式 (6) 中 第 ;个 乘积 UE ÀG.s) Netw 中 有 -2 个 中 间 节 点 与 et 中 的 中 间 节 点 完全 重合 ， 


表示 第 ?个 乘积 项 中 变量 名 (2) 的 次 数 ， 有 日 gis 
2.2 ”模型 学 习 算法 
2.2.1 3E 96 H6 RE 


MTN 参数 可 以 通过 采样 学 习 与 适当 的 学 习 算 法 来 获得 。 
u(k)eR' 为 系统 输入 ， 


HRR T - LK). y(c-D) ja ， 
x(k*DeR 为 系统 输出 。 误 差 变 化 方程 如 下 : 


i-a xo *D-J(k4D) 


-yXo« 十 了 一 2 sofle" n 


qo) 
设 CcC(k)= (II g^, [Tif K), Jre) 


Â = (6(0),6(),..,6(N -1)) 


中 间 节 点 相等 。 


根据 引 理 3， 当 训练 样 


本 充分 多 时 ， 


而 余下 的 两 个 重合 的 中 间 节 点 县 加 后 与 Ne 中 最 后 一 个 


中 间 节 点 的 衰减 


和 重合 是 必然 的 。 因 此 , 有 必要 使 用 前 


枝 算法 来 精简 网 络 ， 


在 去 掉 风 和 


p 


"n 


网 络 结构 。 


有 剪 枝 特性 Po30， 


节点 后 的 网 络 仍然 可 以 表示 一 个 目标 函数 。 
本 文通 过 正则 化 实现 对 MTN IERI AA 


除 来 精简 MTN 


姑 此 一 些 元 余 的 输 


$5 于 该 正则 化 项 
出 权 值 将 逐步 衰减 到 零 
基于 


;从 而 可 以 删除 与 


之 相连 的 中 间 节 点 ， 


此 ， 本 节 在 CG 算法 的 基础 上 引入 带 有 高 斯 正则 化 项 的 


BUE TERR TE Jy MTN 辨识 模型 的 学 习 算法 。 


标 函 数 : 


JO) = E(wi)+yC(wi) 


本 文采 取 如 下 目 


(12) 


Ep, Ew 通常 为 网 络 误差 平方 和 ; 


2w? Iw 


等 式 右 边 第 二 项 为 正则 


9,00 = (aR jm 00) ELR. 化 项 ; 了 为 正则 化 系数 , H720; cuu H wo IEE 
^ ^ ^ ^ 值 。 当 某 些 中 间 节 点 的 输出 权 值 满足 式 (13) 时 , 去 掉 中 间 节 点 。 
Ww, —(W(k),w.(K)...., Ww; (k)) i w| € Wain í (13) 
则 ， 式 (7) 可 以 写成 如 下 形式 ; 其 中 ，w 为 临界 权 值 。 由 于 正则 化 系数 7 的 取 值 对 学 习 结果 
E. dex 有 很 大 影响 , 因此 提出 一 种 在 学 习 过 程 中 动态 改变 7 的 方案 ， 
ui ai E 即 在 每 次 数据 中 心 修正 后 , 动态 修改 7 。MTN 辨识 模型 w 的 
1 调整 方法 采用 如 下 等 式 ?339。 
N’ poo d (8) wi(K+D=wi (K+A p. +yC(w) (14) 
Der (13112 tw (k^) KB, ms, BOEXODHPÉBSXH: 7, coo 与 式 (12) 中 的 
含义 相同 。 
WE EAT "BM, S KWAR Weigend 等 人 P2330 给 出 的 正则 化 参数 调整 方法 ， 
ab 20 ostii. QN n 在 学 习 过 程 中 随时 检测 以 下 误差 量 之 间 的 关系 : 
£5 2 Ny (Ger), -Y yk+Dek ; 
em 2,8 XT (PW N22 + Dé(k) ©) Ek): 当前 时 刻 权 值 调节 时 的 误差 ; 
& gg. p-2Fewrw), i--2 S yaw EG -D : 前 一 时 刻 权 值 调节 时 的 误差， 
9o o Ov, ' x N 各 : ) ， N 名 ? Alk) : 当 少 条 护 的 加 权 平 淘 误 党 ENHA) 9 A(c- D 0—:) EQ) ， 
则 ， 可 以 得 到 其 中 < 为 接近 于 1 的 系数 ; 
pon (10) D: 期 望 误差 值 ， 如 果 没 有 先 验 知识 ， 可 设 定 D-0, 但 
WE y-G0.0.—»0) , 0, ESL E 可 以 重新 写成 如 下 形式 : ”此 时 计算 时 间 可 能 较 长 。 
WT PES c-r 对 7 进行 动态 调节 ， 具 体 规则 如 下 : 
di d E LI NAA moe 如 果 EGO«EX-D, Bü EQO-D, WjryQO- —DAy ; 
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如 果 ED>EK-D ， 
X(k) 2 y(k-D—-Ay ; 

如 果 EGO0ZE(Kk-D ， 或 EWA) ， 
OO=oxk-D， 其 中 为 接近 1 的 系数 。 
2.3 MTN 辨识 方案 
本 文 针对 MIMO 非 线 性 时 灌 系 统 ， 提 出 基于 MTN 
的 辨识 方案 。 HI, MTN 作为 辨识 模型 ; 其 次 ,WE-CG 
算法 作为 MTN 辨识 模型 的 学 习 算 法 ; WE 算法 对 MTN 
辨识 模型 结构 进行 精简 从 而 降低 计算 复杂 度 以 满足 实 
时 性 要 求 。 具 体 过 程 如 算法 1 所 示 。 

算法 1 基于 MTN 的 MIMO 非 线 性 时 滞 系 统 辨 识 

步骤 1: 确 定 MTN 辨识 模型 结构 与 初始 权 值 ; 
步骤 2: 按 照 误差 标准 函数 (7)，MTN 辨识 模型 
步骤 3: 由 式 (14) 更 新 wi ; 

步 又 4: 重 复 步骤 2-4， 直 到 误差 满足 给 定 精 度 停 止 训练 。 


3 ”实验 结果 及 分 析 


本 文通 过 一 个 数值 仿真 例子 和 一 个 实例 -连续 搅拌 反映 
器 系统 (continuous stirred tank reactor, CSTR) 过 程 ， 来 验证 所 
提 方 案 的 准确 性 与 实时 性 。 
3.1 算 例 

考虑 如 下 的 MIMO 非 线性 时 滞 系 统 : 
2y(k - Du? (k -1) 


或 EQ)«AQ0 , H EGO2D , W 


H ED>Z , Wi 


WE -CG 算法 进行 训练 ; 


[RS F0.3u,(k 1— d) &0.7u (k —2)4- 
Ge D A 10386717) 070-2) 
0.2u, (k — d esin(0.5u, (k —d)) 
e 2 E 
do LRRD a Su (k-1-d)+0.5u(k-2)+ 
12 yik-D 
0.2u, (k — d) + sin(0.Su, (k — d —2)) 
HEH, d-6. 


输入 为 单位 阶 跃 信号 ,MTN 辨识 模型 采用 6-28-2 的 
结构 ， 即 6 个 输入 节点 , 展开 到 2 次 ，WE-CG 算法 作为 
其 学 习 算法 。CG 学 习 算 法 设 定 如 下 : ERRAN 200. 
WE 算法 设 定 如 下 : 参考 权 值 为 0.1， 正 则 化 系数 初 值 为 
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denn 


(a) 辨识 结果 (b) 学 习 误 差 曲 线 


106 
时 间 


(c) 正则 化 系数 学 习 变化 曲线 


(d) 网 络 剪 枝 过 程 
图 3 子 系统 2 辨识 
Fig.3 Identification of sub-system 2 


2 与 3 为 辨识 过 程 ， 其 中 ， 图 2(a) 和 图 3(a) 为 辨识 结 
果 ， 图 2(b) 和 图 3(b) 为 学 习 误 差 曲线 ， 图 2(c) 和 图 3(c) 为 正 
则 化 系数 的 学 习 变 化 曲线 , 图 2(d) 和 图 3(d) 为 网 络 剪 枝 过 程 。 
尽管 正则 化 系数 的 变化 并 不 是 平滑 的 学 习 过 程 ， 但 WE 算法 
的 平滑 的 学 习 曲 线 表 明 MTN 权重 变化 过 程 不 会 导致 严重 波 
动 的 训练 误差， 并 对 最 后 的 剪 枝 的 结果 影响 较 小 。 

从 仿真 结果 中 本 文 可 以 看 到 ， 针 对 第 一 个 子 系统 ，WE- 
CG-MTN 方案 和 CG-MTN 方案 训练 的 均 方 误差 分 别 为 
0.0318 和 0.0285， 测 试 的 均 方 误差 分 别 为 0.0365 和 0.0364; 
对 于 第 二 个 子 系统 , WE-CG-MTN 方案 和 CG-MTN 方案 训练 
的 均 方 误差 分 别 为 0.0172 和 0.0172, 测试 的 均 方 误差 分 别 为 
0.0521 和 0.0520。 可 以 看 出 ，WE 算法 的 引入 虽然 对 学 习 误 
差 有 少许 影响 , 但 是 对 于 两 个 子 系统 ,去 除了 元 余 中 间 节 点 后 ， 


0， 滤 波 系数 为 心 0092 ，y=0.995， 正 则 化 系数 的 增 量 为 
0.5x10* , 1 ÉRIC AR BUE I BIB 0.05. HRH 200 个 采 
样 时 刻 进 行 训 练 ，100 个 采样 时 刻 用 于 测试 。 

在 图 2、3 中 , yl 和 y2 为 系统 输出 ，WE-CG-MTN 为 所 
提 辨 识 方案 ，CG-MTN 为 对 比方 案 。 


m 


30 — 100 — 120 40 160 18) — 200 2 w 60 w 10 Ho 1 
时 间 人) 时 间 00 


(a) 辨识 结果 


(c) 正则 化 系数 学 习 变 化 曲线 (d) 网 络 剪 枝 过 程 
图 2 子 系统 1 辨识 
Fig.2 Identification of sub-system 1 


100 
ato 


其 节点 分 别 只 有 11 个 和 19 个， 这 大 大 减少 了 计算 量 ， 增 加 
了 实时 性 能 .WE-CG-MTN 5 CG-MTN 辨识 方案 的 复杂 度 分 
析 如 表 2 所 示 。 


表 2 复杂 度 分 析 


Tab.2 Computation complexity analysis 


预测 模型 构建 方案 结构 ”加 法 乘法 浮 点 数 总 和 
WE-CG-MTN 61122 20 22 42 
CG-MTN 6-28-2 54 56 110 


从 表 2 看 到 经 过 精简 后 的 MTN 辨识 模型 的 浮 点 数 远 远 
小 于 经 典 的 MTN 辨识 方案 ， 也 就 是 说 本 文 的 计算 量 远 远 小 
于 经 典 的 MTN 辨识 方案 , 即 所 提 方 案 具有 很 好 的 实时 性 能 。 

需要 指出 的 是 , MTN 相当 于 神经 网 络 的 一 个 神经 元 , 其 
计算 量 远 远 小 于 神经 网 络 ， 相 关 文 献 已 经 对 传统 MTN 辨识 
模型 和 BP 神经 网 络 辨识 模型 的 准确 性 与 实时 性 进行 了 对 比 
分 析 ， 表 明 传统 MTN 辨识 模型 优 于 BP 神经 网 络 辨识 模型 
09]。 这 里 ， 本 文 对 比 了 所 提 方 案 与 传统 MTN 辨识 方案 的 准 
确 性 与 实时 性 分 析 , 通过 图 2、3 与 表 2 看 出 ， 从 准确 性 分 析 
角度 ， 无 论 是 训练 的 均 方 误差 ， 还 是 测试 的 均 方 误差 ， 本 文 
所 提 辨 识 方案 与 传统 MTN 辨识 方案 几乎 属于 同一 数量 级 ， 
而 从 计算 复杂 度 角 度 ， 本 文 所 提 辨 识 方案 的 计算 量 少 于 传统 
的 MTN 辨识 方案 ， 即 实时 性 更 优 。 
3.2 ”实例 -CSTR 过 程 

CSTR 为 大 氧 处 理 过 程 ， 其 为 一 个 由 原油 与 微生物 完全 
混合 的 状态 B139， 如 图 4 所 示 。CSTR 生成 装置 是 带 有 搅拌 
装置 的 密闭 水 箱 , 其 是 一 个 具有 时 滞 的 典型 非 线 性 化 工 过 程 。 


m 


T 
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CSTR 模型 采用 如 下 模型 1; 


da _ - (t) 
a AOD Ole (Je 


Clu (t7) - x) 


)+ 


x(t) 


=x) + BD, U- x (Dlexp( oe (Jo 


Clu (t - 7) 7x40] 
离散 化 ， 得 
x(O=xk-D-xKk-DAT+DU_xK-D] 
wn(k-D) 
1+x (k -D/o 


exp( JAT + C[u (k - d) - x, ()]AT 


x(k) 2 x(k -1) x, (k -DAT + BD,[1— x, (k — 1] 
nk-l) 
I*x,(k-D/o 


rH, D,-0008, o-10, B-6, C-05, d-6, 采样 时 间 AT 


exp( JAT + C[u, (k — d) — x, (]AT 


为 0.1s。 


结 


PK)=AK); 


»((k)-x(K) , 


物料 


图 4 CSTR 过 程 
Fig.4 CSTR process 
输入 采用 单位 阶 跃 信 号 ，MTN 辨识 模型 采用 6-28-2 的 
构 ， 即 6 个 输入 节点 , 展开 到 2 次 , WE-CG 算法 作为 其 学 


习 算 法 。CG 学 习 算 法 设 定 如 下 : 迭代 次 数 为 1000. WE 算 


结果 ， 图 5(b) 和 图 
正则 化 系数 的 学 习 变 化 


算法 的 平滑 的 学 习 | 


法 设 定 如 下 : 参考 权 值 为 0.1， 正 则 化 系数 初 值 为 0， 滤 波 系 
数 为 x=0.092 , 
元 余 权 重 的 阔 值 为 0.05。 其 中 用 1000 个 采样 时 刻 进行 训练 ， 
200 个 采样 时 刻 用 于 测试 。 


7=0.995 ， 正 则 化 系数 的 增 量 为 05x10* ， 消 除 


图 5 与 图 


[91 


为 辨识 过 程 ， 其 中 ， 图 5(a) 和 图 6(a) 为 辨识 
6(b) 为 学 习 误 差 曲 线 ， 图 5(c) 和 图 6(c) 为 
线 ， 图 5(d) 和 图 6(d) 为 网 络 剪 枝 过 
。 尽 管 正则 化 系数 的 变化 并 不 是 平滑 的 学 习 过 程 ， 但 WE 
线 表 明 MTN 权重 变化 过 程 不 会 导致 严 
并 对 最 后 的 剪 枝 的 结果 影响 较 小 。 

从 仿真 结果 中 本 文 可 以 看 到 ， 针 对 第 一 个 子 系统 ，WE- 


CG-MTN 方案 训练 的 均 方 误差 为 0.0024, 测试 的 均 方 误差 为 
0.0052， 对 于 第 二 个 子 系统 ，WE-CG-MTN 方案 训练 的 均 方 
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间 层 节点 数 为 28 个 。 当 输入 节点 数 不 变 , 随 着 展开 次 数 的 增 
加 ,中间 节点 数 会 随 之 增加 , 造成 节点 见 余 的 发 生 。 对 于 WE 


参数 的 选取 ， 其 决定 了 剪 枝 的 效果 ， 通 常 通过 多 次 实验 来 确 
定 参数 的 选择 。 本 文 看 到 在 引入 WE 算法 后 ， 去 除了 中 间 元 
余 节点 ， 减 少 了 计算 量 。 但 是 1 
的 时 间 来 做 实时 性 分 析 是 不 可 能 的 ， 因 此 本 文 利 ) 
估计 方案 的 计算 复杂 度 ， 


Z 


于 使 用 相同 的 计算 机 ， 相 同 
个 乘 号 都 叫 作 一 个 浮 


个 加 号 和 


点 ， 加 号 和 乘 号 的 总 数 作 为 总 的 浮 点 数 ， 即 浮 点 数 的 总 和 代 
表 着 计算 复杂 度 。 本 文 对 MTN 模型 (CG-MTN) 和 精简 MTN 


模型 (WE-C 


G-MTN) 的 浮 点 数 进行 了 对 比 , 说 明 精 简 MTN 模 


型 浮 点 数 远 远 小 于 MTN 模型 。 同 时 在 以 往 研究 中 分 析 了 
MTN 模型 浮 点 数 远 远 小 于 神经 网 络 模 型 。 这 些 说 明 本 文 所 提 
较 好 的 实时 性 能 。 


方案 有 较 低 的 计算 复杂 度 ， 拥 有 


Jf 


(a) 辨识 结果 


(c) 正则 化 系数 学 习 变化 曲线 


70 


3e) S00 od 
EIS 


(d) 网 络 剪 村 过程 
子 系统 1 辨识 
Identification of sub-system 1 


图 5 


Fig. 5 


1 26 300 400 0500 
时 间 w 


(c) 正则 化 系数 学 习 变化 曲线 


ETE] 
Mgr ó 


(d) 网 络 剪 枝 过 程 
子 系统 2 辨识 


图 6 


误差 为 0.0065， 测 试 的 均 方 误差 分 别 为 0. 


0108。 可 以 看 出 ， 


WE 算法 的 引入 虽然 对 学 习 误 差 有 少许 影 


个 


H, E&x F 


子 系统 ， 去 除了 元 余 中 间 节 点 后 ,其 节点 分 别 只 


2 个 和 7 个 ， 


这 大 大 减少 了 计算 量 ， 增 加 了 实时 性 能 。 


需要 说 明 的 是 MTN 辨识 模型 结构 的 选择 ， 本 文采 用 6- 


28-2 的 结构 ， 即 6 个 输入 节点 ， 展 开 到 2 次 ， 由 图 1 可 知 中 


Fig. 6 
4 ”结束 语 


针对 MIMO 非 线性 时 滞 系 统 提出 了 基于 MTN 的 辨识 方案 。 
a) 利用 MTN 的 强大 逼近 性 ，MTN 作为 辨识 模型 。 
b) WE-CG 算法 作为 MTN 辨识 模型 的 学 习 算法 ，WE 算 
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法 可 以 有 效 精简 MTN 辨识 模型 结构 ,从 而 降低 算法 复杂 度 ， 
提高 模型 实时 性 能 。 

实验 结果 表明 , 本 文 所 提 辨 识 方案 够 准确 地 对 MIMO 非 
线性 时 湾 系统 进行 辨识 。 同时, 相 比 传统 的 MTN 辨识 方案 ， 
所 提 辨 识 方案 结构 更 精简 ， 具 有 更 低 的 算法 复杂 度 ， 实 时 性 
更 好 ， 为 工程 实践 和 硬件 实施 提供 方便 。 
虽然 所 提 辨 识 方法 展现 出 良好 的 性 能 ， 但 仍 需 进一步 的 
WEG. 比如: a) 研 究 其 他 的 精简 辨识 方案 ; b) 研 究 多 入 多 出 非 
线性 时 清 系 统 控制 问题 ;ec) 研 究 该 方案 应 用 到 其 他 工程 领域 


ead 


问题 。 
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